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We examine synchronization between identical chaotic sys¬ 
tems. A rigorous criteria is presented which, if satisfied, guar¬ 
antees that the coupling produces linearly stable synchronous 
motion. The criteria can also be used to design couplings that 
lead to stable synchronous motion. Analytical results from a 
dynamical system are presented. 

05.45.+b 

Synchronization between chaotic systems has been the 
subject of many theoretical papers over the last few years. 
It has also been experimentally observed in many sys¬ 
tems ti- Despite these efforts many key issues re¬ 
main open, and there are few rigorous results that en¬ 
sure the stability of synchronous chaotic motion. In most 
cases rigorous results are obtained using Lyapunov func¬ 
tions (!(] [u| . Unfortunately, this method is not a regular 
approach and, in practice, it can only be applied to par¬ 
ticular examples. The other rigorous approach is that 
of Ashwin et. al. JlJ]. To apply this approach one must 
show that all normal Lyapunov exponents are negative 
for all measures of the dynamics. For typical dynamical 
systems this leads to an intensive numerical analysis. 

There are a few special types of coupling between non¬ 
linear systems where rigirous analysis of the stability of 
synchronization is straightforward. One type is when the 
coupling transforms the driven system into a stable lin¬ 
ear system with time dependent driving. A second is 
when the coupling is diagonal between all of the vari¬ 
ables flj]. In many practical cases these types of coupling 
can’t be achieved. Thus, the state of the art does not give 
a practical answer to the following important question: 
Given an arbitrary dynamical system how can one design 
a physically available coupling scheme that is guaranteed 
to produce stable synchronous chaotic motion? 

This paper examines synchronization between identi¬ 
cal systems with drive/response coupling. The major 
result is a rigorous criteria which, if satisfied, guaran¬ 
tees linearly stable synchronous motion. More impor¬ 
tantly, the criteria can be used to design couplings that 
produce stable synchronized behavior. The criteria only 
uses knowledge of the uncoupled dynamics, and many 
of the important calculations can be performed analyti¬ 
cally. Furthermore, the linearized stability equations we 
examine arise in many other problems that have recently 
appeared in the literature. A discussion of this last issue 


is in our longer manuscript |]l5|| . 

Drive response synchronization between identical sys¬ 
tems is modeled by 


F(x;f) 

(1) 

F(y; t) + E(x — y), 

(2) 


where x is driving dynamics, y is the response dynam¬ 
ics, and E is a vector function representing the coupling. 
For these equations x, y G lR d and E(0) = 0. Synchro¬ 
nization occurs on an invariant manifold given by x = y. 
Obviously, if the coupling strength is below some critical 
threshold then stable synchronous motion will not occur. 
For some (F,E) pairs stable synchronous motion occurs 
only within a finite range of coupling strengths while for 
others synchronization is never stable. 

If one defines deviations from synchronization by w e 
y — x then Eqs. © and (||) lead to the following lin¬ 
earized equation for motion transverse to the synchro¬ 
nization manifold 

^ = [DF(x)-DE(0)]w. (3) 

In this equation DF(x) is the Jacobian of F evaluated 
on the driving trajectory, x, and DE(0) is the Jacobian 
of E evaluated at 0. The synchronization manifold is lin¬ 
early stable if lim^oo ||w(f)|| = 0 for all possible driving 
trajectories x(f) associated with the chaotic attractor of 
the driving system. 

To determine the behavior of w (t) in this limit divide 
DF(x) — DE(0) into a time independent part, A, and 
an explicitly time dependent part, B, 

DF(x) - DE(0) = A + B(x). 

(The nonuniqueness of this decomposition will be re¬ 
solved later.) Assume A can be diagonalized, transform 
to the coordinate system defined by the eigenvectors of 
A, and rewrite the linearized equations of motion as the 
following integral equation |l5| 

z(f) = U(f, to)z(t 0 ) 

+ 1 U(f, s) [P _1 B[x(s)]P] z(s)ds. (4) 

Jt 0 

In this equation z = P _1 w where P = [ei e 2 ■ • ■ e^] and 
ei, e 2 , ■ ■ - ed are the eigenvectors of A. The ordering of 


the eigenvectors is given by the corresponding eigenval¬ 
ues, 9?[Ai] > 3i[A 2 ] > ■•• > 3?[Ad], where 3?[A] is the 
real part of A. Also, JJ(t,to) = exp[D(t — f 0 )] is a time 
evolution operator, where D = P : AP is diagonal by 
assumption. 

Linear stability of the synchronization manifold is de¬ 
termined by ||z(t)|| in the t —■> oo limit. If one uses norms 
to convert Eq. (||) into an inequality, and applies Gron- 
wall’s theorem, then one can define the following decom¬ 
position 


A = (DF) - DE (5) 

B = DF (DF), (6) 

where (•) denotes a time average along the driving tra¬ 
jectory. In terms of this decomposition the criteria for 
linear stability of synchronous motion is j^] 

- K[Ar] > (HP’ 1 |DF(x) - (DF)] P||) . (7) 

Equations (||) (|7;) are our major results. They repre¬ 
sent definitions and conditions that indicate when syn¬ 
chronous motion along a particular driving trajectory is 
guaranteed to be stable to small perturbations in direc¬ 
tions transverse to the synchronization manifold. The 
criterion is rigorous and sufficient. However, because it 
is based on norms it is not necessary. Indeed, numer¬ 
ical experiments indicate that it tends to overestimate 
the necessary coupling strengths (l5|. Also, since the 
integral in Eq. (|7i) is positive semi-definite the inequal¬ 
ity can’t be satisfied unless 3?[Ai] < 0. This condition 
is reminiscent of the discussion of conditional Lyapunov 
exponents found in previous references. 

The decomposition in Eqs. (|j) and (||) is optimal in 
the sense that it minimizes the right hand side of Eq. (|7j). 
We speculate that minimizing this integral gives one the 
chance at satisfying the inequality. Furthermore, by in¬ 
serting Eq. (i into a Volterra expansion of Eq. (® one 
can show that, to second order, the criteria for linear sta¬ 
bility is 3FJ[A]_] < 0 [|n|. For any other decomposition this 
approximate stability criteria will be correct to only first 
order. Finally, for this decomposition Eq. ([?]) reduces to 
3?[Ai] < 0 for fixed points, a result that will not hold for 
other decompositions. 

Equations (®-( 0 ) depend explicitly on the measure of 
the driving trajectory. Gupte and Amaritkar examined 
synchronization using unstable periodic orbits as driv¬ 
ing trajectories |J|. This, and later papers, show that, 
for fixed coupling strength, different driving trajectories 
have different stabilities Q,[L5],[F7|. Recently, Hunt and 
Ott ]l8} ] numerically examined time averages on different 
measures of a chaotic dynamical system and found that 
they tend to assume their largest values on unstable pe¬ 
riodic orbits with the shortest periods. This behavior is 
also discussed in Ref. Ill- 

Given these observations we conjecture that, in many 
practical cases, the unstable fixed points of the driving 
system will be the first measures on the synchronization 


manifold to go linearly unstable as the coupling strength 
is changed. Thus, these trajectories should be the first 
measure to check for linear instability. This conjecture is 
examined in our longer paper and is found to be true for 
the examples studied p~5| . 

Equation ( 0 ) has a geometrical interpretation which 
can be used to design couplings that yield stable syn¬ 
chronous motion. The elements of DE(0) define a pa¬ 
rameter space and each side of Eq. (j7|) defines a function 
in this parameter space. Thus, £ x and £a, respectively 
defined by (HP- 1 [DF(x) - (DF)]P||) = const. = C\ 
and — 3f[Ai] = const. = C 2 , are families of surfaces in 
this parameter space. The boundary of the portion of 
the parameter space that yields linearly stable synchro¬ 
nization is the intersection of these families of surfaces. 
By choosing the elements of DE(0) on portions of Ea 
that are “above” £ x one insures that the poles of A are 
sufficiently far into the left half plane to insure stability. 
Thus, designing a coupling is similar to pole placement 
in control theory El- 

As an example we present an analysis of the following 
dynamical system studied by Ott and Sommerer [|l7| 


dx 

dt 

dv x 

dt 

dy_ 

dt 

dv y 

dt 


= v x 

= -vv x + Ax (l - X 2 ) + y 2 + fo sin (cut) 
= 2 Vy 

= —vv y — 2y (x — p) — Aky 3 


( 8 ) 


where v = 0.05, fo = 2.3, w = 3.5, k = 0.0075 and 
p = —1.5. Originally, Ott and Sommerer examined the 
stability of the invariant manifold defined by y = v y = 
0. Their results indicate that for these parameter values 
motion on this manifold is chaotic, the manifold itself is 
unstable, and only one stable attracting set exists in 1R 4 . 

As before, x denotes the driving system and y denotes 
the response system. In principle the driving trajectory, 
x = [x, v x , y, v y ] £ 1R 4 . However, for this example we 
consider a driving trajectory restricted to the manifold 
examined by Ott and Sommerer. For this type of driving 
DF (x) assumes a block diagonal form. If we use block 
diagonal coupling then Eq. (|j) decomposes into motion 
parallel to, and perpendicular to the manifold examined 
by Ott and Sommerer. 

For perpendicular motion 


dwW 

dt 


DF (i) (x) 


DE (i) (0) 




( 9 ) 


where 


DF (1) (x) = 
DE (i) (0) = 
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FIG. 1. Period 1 and period 2 orbits of the Ott-Sommerer 
model. 


and </(-*-)(x) = —2 (x — p). (Equation (j9l) is the same 
linear stability equation studied by Ott and Sommerer, 
Eqs. (7) and (8) in Ref. |0). An equation similar to 
Eq. (||) involving DF^", DE^\ and^^^(x) = 4(1— 
Zx 2 ) exists for motion parallel to the manifold. (For the 
remainder of this letter we drop the _L and || superscripts, 
and trust the reader to perform calculations in both the 
perpendicular and parallel subspaces.) 

It is easy to show that the eigenvalues of A are 

= —{v + ei + e 2 ) 

* 2 

± \ [i v + £ 2 - ei) 2 + 4(1 - e 4 )((ff> - e 3 )] 1/2 . (10) 

If A± are complex then —5R[Ai] can be made arbitrarily 
large by increasing ei and/or £ 2 - The case for real A± 
is more complicated, however, numerical results indicate 
that — Jf[Ai] is maximized when A± are complex [^ 5 ). Be¬ 
cause, (||P -1 [DF(x) — (DF)]P||) diverges as A± tran¬ 
sitions from real to complex e’s associated with this tran¬ 
sition should to be avoided. These observations suggest 
that in order to satisfy the condition for linear stability 
of the synchronization manifold one should choose e’s so 
that A± are complex with imaginary parts that are not 
near zero. 

It is possible to show that if A± are complex then the 
condition for linear stability of synchronous motion is 

v + ei + e 2 > 4C (\Ag\) , (11) 


where 


C = 


{v + £2 


- (1 - e 4 ) 2 _ 

£i) 2 + 4(1 — e 4 )((< 7 ) -£ 3 ) 


1/2 


( 12 ) 


Measure Type 



(lAff^l) 

<|A^|) 

Period 1 

-1.223 

-6.307 

0.4769 

5.142 

Period 2 

-3 

-7.767 

1.678 

10.30 

SBR 

-3 

-7.038 

1.714 

7.856 


TABLE I. 


and A g(x) = g(x ) — (g). Equations ( |l0| ) (p~2|), and the 
conjecture that the e’s should be chosen so that A± are 
complex, is an analytic solution to the rigorous criteria 
for synchronization. Equating C to a real positive con¬ 
stant, in effect, selects surfaces from the families E x and 
Ea- Each driving trajectory, x, corresponds to a different 
surface. 

Since Eqs. (||) do not have fixed points we examined the 
SBR measure, and the measures associated with the peri¬ 
odic orbits shown in Fig. 0 (the SBR measure is shown in 
Fig. 1 of Ref Jn|). Table g shows numerically calculated 
values for {g) and (|Agr|). 

We now explicitly examine several types of driving. 
The first is diagonal driving. It is the first of the spe¬ 
cial cases where rigirous results are straightforward |lj,g . 
Diagonal driving uses all components of x and chooses 
£3 = £4 = 0, £i = £2 = £ 0^. The parameter space is 
the real line, IR. For this type of driving A± are com¬ 
plex on each of our measures for all values of e. Since 
C is independent of £ its value is fixed, E x is the en¬ 
tire parameter space, and Ea is a family of points in IR. 
Thus, for a particular driving trajectory the boundary 
for linear stability of the synchronization manifold (the 
intersection of Ea with Ex) is given by a point in IR. 
The rigorous criteria for synchronization, Eq. ©, is 


£ > ~ 2 ^ 


-1 


1/2 


JZ 2 + 4 (g) 


Driving via postion uses only the position variables, x 
and y. The simplest example is £2 = £3 = £4 = 0 and the 
parameter space is again IR. It is useful to define new 
parameters u = £i + v and w = 1/C. In terms of u and 
w Eqs. ( 0 ) and (jl^) are 


uw > 4 (|Ag|) 

—4 (g) = (u — 2v) 2 + w 2 . 


If (g) < 0 then these equations define a hyperbola and a 
circle, respectively. It is straightforward to show that the 
circle does not intersect the hyperbola on the measures 
we have examined. Thus, for these driving trajectories 
the rigorous condition for synchronization can not be sat¬ 
isfied. (This does not mean that stable synchronization 
will not result from this type of driving. It only means 
that our analysis can not guarantee that stable synchro¬ 
nization will result from this type of driving [l5|] .) 

Another example of this type of driving uses the po¬ 
sitions to drive both the position and the velocity equa¬ 
tions. (This can sometimes synchronize systems when 















simple driving via position does not produce synchro¬ 
nization |^(J.) For this type of driving, e -2 = £4 = 0 and 
the parameter space is JR 2 . Also, A± are complex for all 
C > 0, while Eqs. ( 0 ) and © become 


d > -v + 4C (\Ag\) 
e 3 = i(V-ei) 2 + ( g} + 


4C 2 


These equations define a line and a parabola in JR 2 . 
For any driving trajectory the line and the parabola are 
guaranteed to intersect. Therefore, synchronization to 
any trajectory, x, is guaranteed to be linearly stable 
for coupling strengths on the parabola whose ei value 
is larger than the one associated with the intersection. 

Driving via velocity uses only the velocity variables, 
v x and v y . As an example, let the velocities drive both 
the position and the velocity equations. Thus, £1 = 63 = 
0, the parameter space is JR 2 , and A± are complex for 
C > 0. If we define new parameters u = v + e -2 and 
w = (e 4 — 1)/C then Eqs. ( 0 ) and 0 ) become 

u>4C(\Ag\) 

(2C(g)) 2 =u 2 + (2C(g)-w) 2 . 


These equations define a line and a circle, respectively. It 
is straightforward to show that the circle does not inter¬ 
sect the line on the measures we have examined. There¬ 
fore, the rigorous condition for synchronization can not 
be satisfied on these orbits. 

In this paper we investigated the linear stability of the 
invariant manifold associated with synchronous behavior 
between coupled chaotic systems. Although we explicitly 
examined unidirectional coupling our results are valid for 
bidirectional coupling, and for determining the linear sta¬ 
bility of invariant manifolds within a chaotic system [ p~5| . 

Our major result is the rigorous criteria of Eqs. (^]) . 

When they are satisfied linear stability of synchronous 
motion is guaranteed. The criteria depends on the mea¬ 
sure of the driving dynamics and can yield different re¬ 
sults for different driving trajectories. The criteria can 
also be used to design couplings that produce synchro¬ 
nization between coupled systems. 

In closing this letter we discuss how noise and non¬ 
linear effects influence our results. Assume the driving 
trajectory is a fixed point, x*, and that Eq. ( 0 ) is satis¬ 
fied for e > e* (Equivalently, a period 1 orbit evaluated 
on a surface of section.) For this case the fixed point un¬ 
dergoes a co-dimension one bifurcation (either pitchfork 
or transcritical) at e = e*. Linear stability analysis does 
not take into account the unstable trajectories near x* 
when e > e*. For arbitrarily small noise amplitude there 
exists a range of e values near e* where the noise will 
eventually push the response system beyond one of the 
unstable orbits. When this occurs the response system 
is forced to seek out an attracting state away from the 
synchronization manifold. Also, nonlinear effects could 
cause an unstable orbit to approach x* for some e far 


from e*. If this occurs then small noise levels can also 
result in a loss of synchronization. 
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